Deep Spin-Glass Hysteresis Area Collapse and Scaling in the d = 3 ± J Ising Model 
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We investigate the dissipative loss in the ±J Ising spin glass in three dimensions through the 
scaling of the hysteresis area, for a maximum magnetic field that is equal to the saturation field. 
We perform a systematic analysis for the whole range of the bond randomness as a function of 
the sweep rate, by means of frustration-preserving hard-spin mean field theory. Data collapse 
within the entirety of the spin-glass phase driven adiabatically (i.e., infinitely-slow field variation) 
is found, revealing a power-law scaling of the hysteresis area as a function of the antiferromagnetic 
bond fraction and the temperature. Two dynamic regimes separated by a threshold frequency w c 
characterize the dependence on the sweep rate of the oscillating field. For u) < oj c , the hysteresis area 
is equal to its value in the adiabatic limit uj = 0, while for ui > ui c it increases with the frequency 
through another randomness-dependent power law. 
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Hysteresis in magnetic materials has been a subject 
of interest for quite some time due to applications in 
magnetic memory devices and as a testground for the- 
ories of nonequilibrium phenomena [lHj]. The hystere- 
sis area which measures the magnetic energy loss in the 
material is connected with the Barkhausen noise p| [f| 
due to irreversible avalanche dynamics 0-[I3- Exist- 
ing literature on hysteresis in random magnets focuses 
mostly on random-field models [l2l - fl5| while numerical 
studies on random-bond models are mostly at zero tem- 
perature [T6M22I ] . To our knowledge, there has been no 
finite-temperature study of the hysteresis loss, especially 
in the spin-glass phase where large avalanches are ex- 
pected to be severely prohibited. We here investigate the 
adiabatic and dynamic hysteresis in the the ± J random- 
bond Ising spin glass [23| on a finite, 3-dimensional sim- 
ple cubic lattice with periodic boundary conditions. We 
show that the hysteresis area obeys a scaling relation 
in the whole spin-glass phase, in accord with earlier 
theoretical studies which observe scale invariance over 
the whole range about the critical disorder for various 
disorder-driven systems [l5l - [T7| . Moreover, this scaling 
data collapse is also observed for experimental systems 
over wide ranges of temparature/magnetic- field scaling 
variables: Gingras et al. observed a universal data col- 
lapse for a temparature/magnetic- field scaling variable 
over four decades in a g eometrically frustrated antifer- 
romegnet Y2M02O7 [24]], while Gunnarsson et al. ob- 
served such a data collapse for the short-range Ising spin 
glass Feo.sMno.sTiOa |1- 

The ± J Ising spin- glass model is defined by the dimen- 
sionless Hamiltonian 



(<J> 



(i) 



where (3 = -^Kp is the inverse temperature. The first sum 
in Eq.([T]) is over the pairs of nearest-neighbor sites 
where Jij is the quenched-random local interaction be- 
tween the classical Ising spins s, = ±1. The probability 



distribution function for Jy is given by 

P(Jy) = P S(J lJ + J) + (1 - p)5{J lJ - J) 



(2) 



H in the second term in Eq.([T]) is the uniform ex- 
ternal magnetic field. With a proper choice of units, 
the temperature for the system may be defined as 
T = 1/J. A random distribution of ferromagnetic and 
antiferromagnetic bonds gives rise to frustration and 
yields a spin-glass phase for a range of p values. Ising 
spin-glass models are widely used as a tool for under- 
standing properties of experimental spin glasses such as 
Pro.6Cao.4Mno.96Gafl.n4O3 EH, Fe . 5 Mn . 5 TiO 3 gUll, 
LiHoo.i67Yo.83 3 F4 :2§3, and Cua-zAlMn* |2jj. With- 
out loss of generality we set p < 0.5, since the par- 
tition function is invariant under the transformation 

P,Uf},{sf} "> (1 -P),{sf},{sf}, where A and B 
signify the two sublattices. 

For small values of p and H = 0, the orientational (up- 
down) symmetry is spontaneously broken below a critical 
temperature T c (p) and long-range ferromagnetic order 
sets in. This phase is well understood within the Lan- 
dau picture where the free energy landscape is described 
by two minima at magnetizations ±m(T,p). Beyond a 
critical fraction p c of antiferromagnetic bonds, reducing 
temperature drives the system into a glassy phase. The 
low-temperature phase now retains its orientational sym- 
metry and a new, randomness-dominated phase which 
has a broken replica symmetry appears (30T l3l| . In this 
phase, the free energy landscape is rough, with many 
local minima at significantly nonoverlapping configura- 
tions. Meanwhile, the dynamics slows down to the extent 
that the relaxation time diverges (32|. At high temper- 
atures, T > T c (p), both ordered phases give way to a 
paramagnetic state where the entropic contribution to 
the free energy is dominant. While the critical temper- 
ature strongly depends on p along the ferromagnet-to- 
paramagnet phase boundary, only a weak dependence of 
T c on p is observed for the spin-glass phase [32|, [33| . In 
this study, we investigate the hysteretic behavior of a spin 
glass under uniform magnetic field H that is swept at a 
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FIG. 1. (Color) Hysteresis curves (upper) and relaxation 
times (lower) at high (T = 4.00, left) and low (T = 2.00, right) 
temperatures. Data are for p values either deep in the ferro- 
magnetic (p = 0.05), spin-glass (T = 2.00, p = 0.35) or para- 
magnetic (T = 4.00, p = 0.35) phases, or close to the phase 
boundaries for the ferromagnetic/paramagnetic (T = 4.00, 
p = 0.15) or ferromagnetic/spin-glass (T = 2.00, p = 0.22) 
transitions. For each case, an overlay of 20 distinct runs with 
different random-bond arrangements is shown. 



boundary conditions. We have checked in this study and 
in a previous study [34| that our hard-spin mean-field 
theory results are independent of size for anLxixi 
system for L > 15. A particular realization at a given 
{T,p) is generated by the assignment of the quenched- 
random coupling constants Jy according to the proba- 
bility distribution of Eq. © and, initially, a random and 
unbiased choice of spins Sj = ±1. In order to determine 
the hysteresis curves, the system is first saturated by a 
sufficiently large external field H s , the minimum value 
of H for which Eq.© yields an average magnetization 
m = {^/L 3 )J2i m i — 1 within an accuracy e m = 10~ 6 . 
Then, the path H s — > — H s — > H s is traversed with steps 
AH = H s /100 or smaller. For each incremental change 
of the field, the system is allowed to relax a number of 
time steps r = 1/uj. A time step corresponds to succes- 
sive iterations of Eq.(|3]) on L 3 arbitrarily chosen sites. 
An infinitely slow sweep is obtained as the limit r — > tu, 
where the HSMF equations converge to a self-consistent 
solution within the tolerance interval e m . Thus, tn is the 
relaxation time of the system. 

The infinitely-slow-sweep hysteresis curves obtained in 
the ferromagnetic and spin-glass phases are shown in 
FigjT] The usual jump in the magnetization at a co- 
ercive field H c , observed for small p, is associated with a 
system-wide avalanche in the ferromagnetic phase. For p 
larger than a critical value p c , this picture is replaced by 
a slanted hysteresis curve and a smaller hysteresis area, 



constant rate ui. A past computational study similar to 
ours [34j considers a time-dependent quenched-random 
magnetic field that is conjugate to the spin-glass order 
parameter. 

We use hard-spin mean-field theory (HSMFT), a self- 
consistent field theoretical approach [3J-[5fJ that pre- 
serves the effects due to frustration (crucial for the spin- 
glass phase) generated by the randomly scattered anti- 
ferromagnetic bonds. HSMFT is defined by the refined 
set of self-consistent equations 




JJP(m j ,s i ) 



tanh ^ Jij s j + H 



for the local magnetization rrii at each site i, whose 
nearest-neighbors are labeled by j . The single-site prob- 
ability distribution is 



P(TOj,Sj) 



1 + mjSj 



(4) 



The local magnetization rrn at site i satisfies — 1 < rrii < 
1. The hard-spin mean- field theory Eq.(3) has been dis- 
cussed in detail in Refs. [33 - l50j . 

HSMFT has been successfully applied to spin 
glasses [13, 5H . In this paper we make use of the method 
to investigate the scaling of the hysteresis area under a 
uniform, time-dependent magnetic field. To this end, 
we consider a 20 x 20 x 20 cubic lattice with periodic 




FIG. 2. (Color) Logarithmic contour plot of the infinitely- 
slow-sweep hysteresis area Ao as a function of antiferromag- 
netic bond probability p and temperature T = 1/J. The 
thick vertical line denotes the phase boundary between the 
ferromagnetic and the spin-glass phases as described in the 
text, while the other thick line bounds the paramagnetic 
phase where the infinitely-slow-sweep hysteresis area is less 
than the precision used in the consistent-field calculations, 
i.e., A < 10" 6 . 
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FIG. 3. (Color) Infinitely-slow-sweep hysteresis area Ao, as a function of antiferromagnetic bond probability p for temperatures 
(indicated in the color legend) 1/T = J = 0.10, 0.11, . . ., 0.50 (left) and as a function of temperature T — 1/J for antiferromag- 
netic bond probabilities (indicated in the color legend) p = 0.00, 0.01, . . ., 0.50 (right). Each curve is a 10th degree polynomial 
fit to the averages over 20 realizations. 



typical of spin-glass materials [3, [Tl|, Hi| . This converse 
hysteretic behavior, associated with the Barkhausen 
noise 5, 6], is a consequence of the power-law distribution 
of avalanches which is well established 
Hl HI HI H! HI for several frustrated systems with 
quenched disorder. The hysteresis area disappears in the 
paramagnetic phase. 

In FigfJl we present the infinitely-slow-sweep hystere- 
sis area globally, for all temperatures and antiferromag- 
netic bond probabilities, on a logarithmic color-contour 
plot. The hysteresis area Ao vanishes in the region 
shown in dark blue, which corresponds to the param- 
agnetic phase, while it is nonzero in the ferromagnetic 
and spin-glass phases, respectively on the left and right 
of the lower half of Figl2J The para/ferro and para/spin- 
glass phase boundaries are easily determined by locating 
the temperature at which Aq vanishes (i.e., falls below 
e m ). A set of p scans for different temperatures and a 
set of temperature scans for various p values are given 
in Fig|3J The low-temperature ferro/spin-glass bound- 
ary is located at p c ~ 0.22 and is calculated as the 
inflection point for the maximum slope of the hystere- 
sis curve as a function of antiferromagnetic bond prob- 
ability 16]. The phase boundaries are consistent with 
the well-known phase diagram for the 3-dimensional ± J 
model [33[ and in fair comparison with the experimental 
temperature-concentration phase diagrams of the vari- 
ous Eu^Si'i—^SySei—ji^ solid (o-R 2 )i~x(p-ii2)x, and AuFe 



teresis area in the form Aq = Aq(p, J), where p 



systems reviewed in [32 1 



We here focus on the scaling form of the hysteresis area 
in the spin-glass phase and show that a unique scaling- 
function governs the whole range of p and J within the 
spin-glass phase. To this end, we first express the hys- 



and J 



j-j c 



are the reduced displacements from phase 



boundaries. We then postulate the multivariate scaling 
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FIG. 4. (Color) Scaling of the hysteresis area in the spin- 
glass phase as a function of reduced antiferromagnetic bond 
concentration p and the reduced bond strength J, for various 
p values as shown in the color legend. The scaling function 
f(x) given by the RHS of Eq.© on which all data points 
collapse is consistent with a same power law within the entire 
spin-glass phase. 
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form 

A (p,j) = \ c Ao(\ a p,\ b j), (5) 
which by letting A = p~ x l a reduces to 

A (p,J)=p- c / a Ml,p- b/a J), (6) 

Defining v = c/a, ji = —b/a, and f(x) = Aq(1,x), we 
obtain 

p v A (p,J) = f(p»j). (7) 

The sought collapse is obtained by the choice of scaling 
exponents /i = 1 and v = 2. Data shown in FigJ3] col- 
lapses onto a single curve shown in FigJU where the LHS 
of Eq.Q is plotted against the argument on the RHS for 
28 evenly spaced values of p above p c . The origin corre- 
sponds to the phase boundary between the spin-glass and 
paramagnetic phases. The log-log plot of the same col- 
lapse shown in the inset of Fig 0] suggests that the scaling 
function has the form f(x) oc x 1 - 72 , yielding a hysteresis 
area A oc p a P with a ~ -0.28 and (3 ~ 1.72. Inter- 
estingly, unlike the case of usual critical phenomena, the 
scale-invariance applies to the entire spin-glass phase and 
not just to the vicinity of the critical phase boundary. 

Having analyzed the limit with infinitely-slow-sweep 
rate, we next consider the dynamic hysteretic response as 
a function of the magnetic field frequency. One can sim- 
ulate the finite oscillation frequency by iterating Eq.Q 
for a predetermined number of steps t, instead of waiting 
until a steady-state is reached. The sweep rate lu = l/t is 
proportional to the frequency of the applied field up to a 
material-dependent spin relaxation time. The hysteresis 
area A(u,p 1 J) deviates from the value at infinitely-slow 
sweep Aq = A(lu = 0,p, J) and increases with increas- 
ing sweep rate lu. This can be understood by observing 
that the laggy response of the magnetization to a time- 
varying field inflates the hysteresis curve along the field 
direction. The typical behavior observed in various ex- 
perimental and theoretical ma gnets ( typically pure mag- 
nets or random-field systems) [52j-|56| is 

A{u,p, J) = A Q +g(p, J)lj", (8) 

where b is the sweep-rate exponent. We investigate 
whether the random-bond Ising spin glass obeys a similar 
scaling relation. 

A typical scan of the hysteresis area as a function of 
lu displays two dynamic regimes, separated by a critical 
sweep rate lu c that depends on p, J, and the system size 
(FigJS]). For a sufficiently slowly varying field lu < lu c , 
the area is pinned at the value Aq. In this regime, the 
avalanches that are triggered by an incremental increase 
in the field decay within a period 1/lu or smaller. For 
faster sweeps (lu > lu c ), the increase in the area follows 
the power law in Eq.©, with a p-dependent exponent 
b. In the ferromagnetic phase with weak disorder, the 
two dynamic regimes are separated by a sharp increase 
in the hysteresis area. This transition gets significantly 
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FIG. 5. (Color) Hysteresis area difference A — Aq versus sweep 
rate ui, for temperatures T = 2.0, 1.0, 0.5 from top to bottom 
and for antiferromagnetic bond fractions p — 0.0, 0.1, . . ., 0.5 
as shown in the color legend. 
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FIG. 6. (Color) Sweep-rate exponent b versus antiferromag- 
netic bond fraction p for temperatures T = 2.0, 1.0, and 0.5. 
The dashed curve depicts the general trend of the sweep-rate 
exponent, while the dotted vertical line marks the phase tran- 
sition from ferromagnetic to spin-glass phase. 



smoother in the spin-glass phase, especially far from the 
ferromagnetic/spin-glass boundary. For larger systems, 
one expects ui c to recede and the power-law behavior to 
dominate. 

Figj6] shows the sweep-rate exponent b calculated as 
a function of antiferromagnetic bond fraction p, at fixed 
temperatures T — 1/J = 2.0, 1.0, and 0.5. The hys- 
teresis area is calculated for the sweep rates cj = 1, 0.5, 
0.3, 0.2, 10 _1 , . . ., 10~ 4 at each p value, after averag- 
ing over 10 realizations. The exponent values are ob- 
tained through fits to the data in the regime uj > uj c 
(typically two decades or more), using the functional 
form of Eq.©. The error bars reflect only the scat- 
ter of the data relative to the fit. In the ferromagnetic 
phase, p < p c , we note that the calculated sweep-rate ex- 
ponents lie in an interval of fairly good agreement with 
the various values obtained previously at p — 0, namely 
b = 2/3 [H-ii] and b = 0.52 ± 0.04 [Hi from mean-field 
theory, b = 0.61 15311 from Glauber dynamics simulations, 
b = 0.495 ± 0.005 [Ei] and b = 0.45 0] from Monte Carlo 



simulations. 

In conclusion, we have considered here the ±J Ising 
model under uniform external field and investigated the 
scaling behavior of the saturation hysteresis area (i.e., 
far from the weak- field limit). We observed that the 
phase diagram can be derived from the hysteresis area 
alone and the ferromagnetic/spin-glass phase boundary 
corresponds to the inflection point wrt bond-randomness 
strength p. When adiabatically driven, the area displays 
a data collapse within the entire spin-glass phase for all 
temperatures and p. The scaling function itself has a 
power-law form and the scale- invariance extends far from 
the phase boundary, deep into the spin-glass phase. 

The dynamical response under a fluctuating external 
field is also interesting. We find that, beyond a threshold 
value ui c , the hysteresis area increases as a function of 
the field-sweep rate uj with a non-universal power-law. 
This behavior is not limited to the vicinity of the phase 
transition. The associated exponent is found to be a 
function the randomness strength p. Moreover, this 
function is independent of temperature. In the limit 
of a pure magnet (p — y 0), we observe good agreement 
with the existing literature, despite the fact that the 
earlier theoretical work applies to a weak driving field, 
while we here consider sweeps across saturation limits. 
Fig |5] suggests that, relative to the ferromagnetic phase, 
the spin glass displays an amplified sensitivity to the 
field-sweep rate, again running in apparent contrast 
with the general wisdom that the hysteretic effects are 
suppressed within a spin glass. In fact, we note that 
the increase in the hysteresis area with u is due to the 
magnet's delayed response to the changing field, and a 
signature of the spin-glass phase is the slowing down of 
precisely such relaxation phenomena. 
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